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PURPOSE:  The  calculation  of  longshore  sediment  transport  (LST)  rates  is  a  key  component  of 
most  coastal  engineering  studies.  While  the  LST  process  is  conceptually  simple,  in  practice  the 
development  of  reliable  rates  is  made  difficult  by  problems  associated  with  collecting  accurate  field 
data,  by  limitations  to  model  predictions,  and  by  substantial  variations  of  the  rates  in  time  and  space. 

This  Coastal  and  Hydraulic  Engineering  Technical  Note  (CHETN)  summarizes  the  state  of 
understanding  of  the  influence  of  grain  size  on  surf  zone  sediment  transport  and  is  a  companion  to 
Smith  et  al.  (2004).  This  CHETN  discusses  details  of  bed-load  and  suspended  load  transport,  and  the 
classical  bed-load  regime  is  shown  to  encompass  two  distinct  modes  of  transport.  Four  LST  models 
with  varying  levels  of  complexity  are  discussed  to  show  how  they  incorporate  the  physics  of  grain 
size  variation  and  its  effect  on  the  transport  rate.  In  addition,  a  relationship  between  the  K  coefficient 
in  the  CERC  formula  ( Coastal  Engineering  Manual  (2002),  Section  III-2-3-a)  and  grain  size  is 
presented.  Finally,  some  inconsistencies  between  theory  and  data  are  discussed  in  the  context  of  the 
interrelationship  between  grain  size  and  beach  slope. 

The  final  conclusion  is  repeated  here.  In  general,  an  increase  in  the  median  grain  size  will  decrease 
LST  rates  in  the  surf  zone.  If  a  simple  exponential  relationship  between  transport  rate  and  grain  size 
(D)  is  needed  (Equation  7:  LST  rate  oc  Dn),  the  most  appropriate  value  for  the  exponent  n  should  be 
of  the  order  of  -1,  as  seen  from  Equation  23.  However,  this  tech  note  argues  that  this  is  clearly  a 
simplistic  view  of  surf  zone  sediment  dynamics.  A  more  realistic  (though  still  highly  simplified) 
approach  would  be  that,  for  fine  grain  sediments  ( D50  on  the  order  of  0. 15  to  0.30  mm),  suspended 
load  transport  should  dominate  and  n  should  be  somewhere  within  the  range  of  -0.5  to  -3.0.  For 
coarse  sands  (D50  around  1 .0  mm),  sheetflow  bed-load  transport  should  dominate  (Figure  6),  and  the 
transport  rate  should  be  nearly  independent  of  grain  size  (n  =  0).  For  large  gravel  and  shingle  (D 50 
>  20  mm),  the  dominant  transport  should  be  in  the  IM  (Initial  Motion)  bed-load  regime,  with  n 
within  the  range  of  -0.5  to  -2.0.  To  state  this  another  way,  the  exponent  n  in  Equation  7  is  itself  a 
function  of  grain  size. 

BACKGROUND:  Though  most  researchers  recognize  the  median  grain  size  as  a  parameter  of  first 
order  importance  in  its  effects  on  the  magnitude  of  sediment  transport,  its  variation  is  not  as  easily 
studied  as  the  variation  of  many  of  the  other  primary  parameters.  In  the  field,  it  is  easy  to  obtain  data 
for  a  range  of  values  for  wave  heights  and  periods,  for  instance,  as  these  vary  constantly  in  time. 
However,  the  grain  size  on  a  beach  typically  shows  no  (or  insignificant)  variations  in  time.  Similar 
problems  can  occur  in  the  laboratory.  Changing  the  wave  height  or  period  in  a  laboratory  study  is 
usually  as  simple  as  reprogramming  a  wave  paddle.  However,  usually,  changing  the  grain  size  in  a 
large  model  is  infrequently  done  because  of  the  large  cost  in  both  time  and  money. 
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Consequently,  the  effects  of  grain  size  on  sediment  transport  rates  is  not  as  well  understood  as  for 
some  other  parameters.  Grain  size  has  occasionally  been  incorporated  into  models  to  make  expres¬ 
sions  dimensionally  correct,  rather  than  based  upon  theory  or  data. 

Relationship  of  Grain  Size  to  Fall  Velocity  and  Shear  Stress.  The  fall  velocity  is  the 
tenninal  velocity  reached  by  a  sediment  particle  as  it  falls  through  an  infinite,  quiescent  fluid,  where 
the  fluid  drag  retarding  the  particle  is  equal  to  the  downward  gravitational  pull.  Many  sediment 
transport  relationships  (particularly  suspended  load  relationships)  contain  a  fall  velocity  parameter 
(ny)  with  or  without  a  separate  grain-size  parameter  (D  or  Dso,  the  median  grain  diameter).  Fall 
velocity  increases  with  grain  diameter.  For  silts  and  clays,  ws  varies  as  D2  and  for  gravels  it  varies  as 
D  \  For  sand-size  particles,  ws  falls  in  a  transition  region  where  it  varies  as  the  grain  size  to  a  power 
between  !4  and  2  ( Coastal  Engineering  Manual  2002,  Section  III- 1).  For  convenience,  this 
relationship  can  be  crudely  approximated  as  a  variation  to  the  first  power  for  sand-size  particles,  but 
to  obtain  accurate  values,  a  formula  such  as  those  of  Hallermeier  (198 1)  or  Ahrens  (2000)  should  be 
used. 

An  important  parameter  in  sediment  transport  is  the  shear  stress  (t),  which  is  the  tangential  force  per 
unit  area  that  a  moving  fluid  exerts  on  the  sediment  bed,  which  causes  transport.  In  a  wave 
dominated  environment,  it  is  generally  related  to  the  horizontal  fluid  velocity  just  above  the  top  of 
the  boundary  layer  ( u )  by  the  expression: 

T  =  /2fwPU2  (1) 

where  p  is  the  fluid  density,  and fw  is  the  friction  factor.  The  stress  on  the  bed  is  strongly  affected  by 
the  nature  of  the  flow  within  the  boundary  layer.  The  boundary  layer  flow  can  be  laminar  or 
turbulent;  it  is  affected  by  the  bed  roughness  and  the  presence  of  bed  forms  such  as  ripples;  and  it  is 
substantially  different  in  steady  versus  oscillatory  flow.  Generally,  an  attempt  is  made  to  represent 
these  effects  as  variations  in  the/,’,  tenn.  In  oscillatory  flow,  the  /,',  tenn  is  usually  inversely  related  to 
the  dimensionless  term  ( a/ks ),  where  a  is  the  orbital  excursion  amplitude  and  ks  is  the  bed  roughness. 
For  flat  beds  with  sand  or  larger  grain  sizes,  ks  is  directly  related  to  the  grain  size.  Several  equations 
have  been  proposed  for  this  relationship  (e.g.,  Jonsson  1966;  Swart  1974;  Kamphuis  1975;  Nielsen 
1979).  For  the  purposes  of  this  discussion,  the  important  result  is  that,  for  sand-  and  gravel-sized 
particles,  an  increase  in  the  sediment  grain  size  will  cause  a  small  increase  in  the  shear  stress  on  the 
bed,  which  will  lead  directly  to  a  slightly  greater  sediment  transport  rate  for  the  same  free-stream 
flow  conditions.  While  this  effect  is  important,  it  can  be  easily  masked  by  other  greater  influences 
that  the  size  of  the  sediment  has  on  the  transport  rate,  as  discussed  in  the  following  paragraphs. 

Minimum  and  Maximum  Beach  Sediment  Sizes.  The  lower  grain-size  limit  for  the  material 
composing  most  beaches  is  in  the  0.10-  to  0.15-mm  range.  This  value  is  largely  controlled  by 
ambient  near-bottom  turbulence  levels.  There  are  two  primary  sources  of  turbulence  in  the  surf  zone: 
breaking  wave  turbulence  and  turbulence  generated  within  the  wave  or  wave/current  bottom 
boundary  layer.  The  turbulence  generated  by  breaking  waves  is  substantially  larger,  but  is  generated 
at  the  top  of  the  water  column  and  its  effect  on  the  bed  is  highly  episodic  on  both  wave  and  storm 
time  scales.  Sediment  picked  up  by  breaking  wave  suspension  events  will  at  least  occasionally  have 
enough  time  to  settle  back  to  the  bed  before  the  next  event.  However,  small  grains  (slits  and  clays, 
with  low  fall  velocities)  are  mobilized  by  the  turbulence  generated  in  the  bottom  boundary  layer  and 
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kept  in  suspension  throughout  a  half  wave  cycle.  The  same  thing  occurs  on  next  half  cycle,  and  the 
material  never  settles  to  the  bed.  Instead,  the  material  is  kept  in  suspension  and  slowly  diffuses  or  is 
transported  into  quiescent  waters  elsewhere,  either  offshore  or  into  back  bays  and  estuaries.  This 
process  is  analogous  to  that  which  supports  wash  loads  in  rivers.  Beaches  with  grain  sizes  smaller 
than  about  0. 10  mm  are  usually  only  found  near  large  river  deltas  (such  as  many  sections  of  the  coast 
of  Louisiana),  where  the  system  is  overwhelmed  by  the  presence  of  fine  sediments  and  coarser 
material  is  unavailable. 

Few  beaches  are  composed  of  material  whose  median  grain  size  is  larger  than  100-300  mm.  The 
maximum  size  of  beach  materials  is  governed  by  the  rate  at  which  large  materials  can  be  delivered  to 
a  beach  by  rivers,  eroding  bluffs,  etc.,  compared  with  the  rate  at  which  these  materials  weather 
(ablate)  once  they  are  on  the  beach.  The  weathering  of  cobble-sized  materials  can  be  quite  rapid  in  a 
geological  sense.  In  comparison,  the  fine  quartz  sands  (0.15  to  0.3  mm)  that  compose  most  beaches 
have  so  little  mass  that  they  are  almost  totally  resistant  to  further  weathering.  For  additional 
discussion,  see  Coastal  Engineering  Manual  (2002),  Section  III- 1 . 

SEDIMENT  TRANSPORT  REGIMES:  The  theory  leading  to  the  development  of  the  CERC 
formula  was  based  upon  the  concept  that  waves  deliver  a  certain  amount  of  energy  to  the  surf  zone, 
and  that  a  constant  proportion  of  this  energy  is  available  to  mobilize  and  transport  sediment 
(Bagnold  1963;  Inman  and  Bagnold  1963;  Bagnold  1966).  Though  remarkable,  this  development 
had  been  found  to  have  limitations.  As  waves  break  and  release  energy  across  the  surf  zone,  the 
proportion  of  this  energy  used  to  transport  sediment  is  not  constant,  but  varies  substantially  through 
surf  zone  parameter  space.  The  understanding  that  surf  zone  sediment  transport  is  not  a  simple 
process  that  can  be  described  by  a  simple  one -term  equation  has  led  to  the  development  of  more 
sophisticated  equations  employing  the  concept  of  sediment  transport  regimes. 

Most  of  the  LST  models  in  use  today  are  either  one-tenn  models  that  calculate  the  total  load  or  they 
have  two  tenns  that  partition  the  LST  into  bed-load  and  suspended  load  components.  However,  there 
are  at  least  five  nearshore  sediment  regimes  that  can  realistically  dominate  the  transport  under 
certain  conditions:  two  bed-load  regimes,  the  suspended  load  regime,  the  ripple  regime,  and  the 
swash  regime.  A  complete  nearshore  transport  relationship  would  include  components  for  the  calcu¬ 
lation  of  transport  in  each  of  these  regimes,  or  include  an  analysis  that  indicated  which  were 
important  for  a  given  set  of  input  parameters  and  then  only  calculate  those. 

The  bed-load  and  suspended  load  regimes  are  discussed  in  detail  in  the  following  paragraphs  with 
emphasis  on  how  the  transport  is  affected  by  grain  size.  Ripple  transport  is  generally  considered  as 
being  the  same  as  suspended  load  transport  over  flat  beds,  except  that  ks  is  dominated  by  ripple 
height  over  the  grain  diameter.  However,  significant  differences  are  recognized,  including  suspen¬ 
sion  occurring  primarily  or  exclusively  at  the  ripple  crest,  sediment  trapping  vortices  of  standard  size 
forming  in  the  lee  of  each  ripple  and  dissipating  each  half  cycle,  and  transport  processes  that  differ  in 
space  along  the  length  of  the  ripple.  Ripple  transport  also  is  visually  quite  distinct  from  suspended 
load  transport  on  flat  beds.  To  the  extent  that  this  complex  process  is  understood  and  can  be 
parameterized,  the  variation  of  longshore  transport  rate  with  grain  size  in  this  regime  is  thought  to  be 
similar  to  that  in  the  suspended  load  regime.  In  both,  the  transport  strongly  decreases  with  increasing 
grain  size. 
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The  swash  zone  is  defined  as  the  area  between  the  minimum  and  maximum  wave  uprush  for  a  given 
tide  stage.  In  the  swash  zone,  transport  occurs  through  the  mechanisms  of  bed-load  and  suspended 
load.  However,  several  factors  make  this  region  distinct.  The  shallow  depths  generally  allow  the 
boundary  layer  to  grow  to  the  water  surface,  making  aspects  of  the  transport  more  similar  to  steady 
flow  in  rivers  than  to  typical  unsteady  transport  by  waves.  Though  generally  a  large  percentage  of 
the  wave  energy  has  dissipated  before  reaching  the  swash,  the  small  depths  in  this  region  mean  that 
the  dissipation  at  the  water  surface  occurs  close  to  the  bed.  This  has  the  potential  to  elevate  near-bed 
turbulence  levels,  compared  to  the  levels  further  offshore  in  the  surf  zone.  The  relatively  high 
transport  rates  frequently  seen  in  the  swash  may  be  a  result  of  this  fact.  Flow  into  and  out  of  the  bed 
(percolation)  is  substantially  more  pronounced  in  the  swash  than  in  deeper  parts  of  the  surf  zone. 
This  both  changes  the  hydrodynamics  of  uprush  and  downrush,  and  it  changes  the  shear  stress  at  the 
bed,  leading  directly  to  changes  in  mobilization  and  transport  of  the  sediment.  While  these  effects 
are  understood  qualitatively,  how  they  alter  the  relationship  between  grain  size  and  transport  rate 
needs  substantial  further  work. 

Bed-load  Transport  Regimes.  Bed-load  is  defined  as  that  portion  of  the  moving  sediment  that 
is  supported  by  intergranular  contact,  either  continuously  through  rolling  or  sliding  motion,  or  inter¬ 
mittently  through  saltation  (Bagnold  1956).  This  is  distinct  from  suspended  load,  which  is  that 
portion  of  the  moving  sediment  that  is  supported  by  the  turbulent  fluid  drag.  These  distinctions  are 
useful  for  model  development  but  make  precise  measurements  of  the  two  quantities  difficult,  as  both 
processes  occur  at  the  same  elevation  above  the  bed,  or  may  be  partially  supporting  the  same  grain  at 
any  given  moment. 

Many  bed-load  to  total  load  sediment  transport  models  incorporate  a  threshold  of  motion  term  to 
deal  with  low  transport  rates.  A  more  physically  realistic  approach  is  to  consider  bed-load  transport 
as  occurring  in  two  regimes,  an  IM  bed-load  regime  and  a  sheetflow  bed-load  regime.  In  sheetflow, 
there  is  a  dense  layer  of  sediment  many  grain  diameters  thick  in  motion  just  above  the  stationary 
(flat)  bed.  (In  this  discussion,  sheetflow  will  refer  only  to  bed-load  transport,  and  will  not  include 
material,  which  may  or  may  not  be  carried  in  suspension  above  it.)  This  layer  rises  and  collapses 
with  each  half  cycle  of  fluid  motion.  For  less  intense  transport  (decreasing  x),  this  moving  sediment 
layer  decreases  in  both  thickness  and  velocity.  This  concept  of  a  moving  sheet  of  material  works 
well  until  the  thickness  of  the  layer  decreases  to  the  order  of  a  grain  diameter.  Below  this  value,  the 
transport  layer  thickness  cannot  continue  to  smoothly  decrease.  Instead,  transport  occurs  by  some  of 
the  grains  in  the  surface  layer  moving  substantial  distances  while  others  are  remaining  stationary. 
This  is  the  IM  bed-load  regime. 

Though  there  are  no  obvious  visual  differences  in  these  regimes  (such  as  distinguishes  ripple  trans¬ 
port  from  suspended  load,  for  instance),  there  are  marked  differences  in  the  behavior  of  the  sedi¬ 
ment.  Several  researchers  have  measured  bed-load  transport  in  the  laboratory  using  a  variety  of 
facilities  and  measurement  techniques.  Figure  1  shows  half  cycle  transport  ratios  when  beds  of 
different  sized  sediments  were  subjected  to  the  same  flow  fields.  This  figure  shows  the  amount  of 
bed-load  transport  for  the  smaller  material  divided  by  the  amount  for  the  larger.  The  x-axis  is  the 
Shield’s  parameter,  0: 
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Figure  1 .  Laboratory  data  for  ratio  of  sediment  transport  occurring  over  a  half  wave  cycle  for  a  smaller 
grain  size  divided  by  a  larger  grain  size,  when  both  are  subjected  to  same  flow  field  (modified 
from  King  1991) 

where  xm  is  the  maximum  orbital  shear  stress  equal  to  Zifw  pum  ,  ps  is  the  sediment  density,  and  um  is 
the  maximum  orbital  water  velocity  above  the  bottom  boundary  layer.  The  shear  stress  was 
calculated  as  suggested  in  Nielsen  ( 1 992)  and  Soulsby  ( 1 997)  with  the  bed  roughness  ks  =  2.5 D,  and 
fw  was  calculated  using  the  formula  of  Swart  (1974).  The  D’s  in  the  Shield’s  parameter  were 
calculated  using  the  larger  grain  diameter  for  each  pair  of  sediments.  Plotted  ratios  are  limited  to 
those  where  both  values  were  measured  by  the  same  researchers  using  the  same  apparatus  and 
measurement  technique,  to  reduce  the  impact  of  experimental  error,  and  where  the  flow  fields  for 
both  tests  (excursion  amplitudes,  periods,  and  maximum  velocities)  varied  by  less  than  3  percent. 
Grain  sizes  for  the  data  shown  in  Figure  1  are  listed  in  Table  1.  Data  also  were  excluded  if  the 
suspension  parameter,  b,  was  less  than  1 .5  for  either  grain  size,  indicating  that  a  substantial  fraction 
of  the  transport  was  occurring  as  suspension.  The  dimensionless  suspension  parameter  is: 


k  u* 


where  ws  is  the  sediment  fall  velocity,  calculated  using  Ahrens  (2000),  k  is  von  Karman’s  constant 


=  0.40,  and  u*  is  the  friction  velocity  = 
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Table  1 

Grain  Sizes  for  Figures  1  and  2 

Author 

#  of  Ratios 

Numerator  Grain  Size 

mm 

Denominator  Grain  Size 

mm 

Manohar  (1955) 

11 

0.28 

1.006 

Kalkanis  (1964) 

15 

1.68,  2.18 

2.18,  2.82 

Abou-Seida  (1965) 

2 

0.30,  1.21 

1.21, 2.61 

Sleath  (1978) 

2 

1.89 

4.24 

Sawamoto  and  Yamashita  (1986) 

5 

0.7 

1.6 

King  (1991) 

32 

0.44 

1.1 

Figure  1  shows  a  distinct  break  in  the  trend  of  the  data  occurring  near  0  ~  0.2.  For  larger  values  of 
the  Shield’s  parameter,  the  transport  is  nearly  independent  of  grain  size.  Below  0  ~  0.2,  there  is  a 
strong  decrease  in  transport  rate  with  increasing  grain  size.  In  a  series  of  papers,  Bagnold  (1954, 
1956,  1963,  1966)  has  shown  that  the  shear  stress,  \m,  necessary  to  transport  a  given  immersed 
weight  of  sediment,  P,  is  governed  by  Coulomb’s  law  of  sliding  friction: 

=  tan  cj)  (4) 

where  cj)  is  the  angle  of  repose.  The  immersed  weight  of  sediment  per  unit  area  one-grain  layer  thick 
is: 


p  =  (p,  -p)gD(l~p)  (5) 

where  p  is  the  porosity  of  the  sediment.  In  terms  of  the  Shield’s  parameter,  a  one-grain  diameter 
thick  layer  corresponds  to: 

0  =  (l-/>)tan(j)  (6) 

Using p  =  0.4  and  tan  cj)  =  0.58,  Equation  6  yields  a  value  of  0  =  0.35  for  a  transport  layer  one  grain 
diameter  thick.  Sawamoto  and  Yamashita  (1987,  their  figure  13a)  present  a  graphical  relationship 
between  the  depth  of  the  moving  sediment  layer  and  the  Shields  parameter.  Their  data  show  that  a 
moving  sediment  thickness  of  a  one  grain  layer  occurs  near  0-0.2. 

Figure  1  indicates  that  in  the  sheetflow  bed-load  regime,  which  occurs  for  0  >  0.2,  the  transport  rate 
is  largely  independent  of  grain  size.  In  this  regime,  the  shear  stress  at  the  bed  is  able  to  mobilize  a 
large  mass  of  sediment  per  unit  surface  area,  and  it  makes  little  difference  whether  that  mass  is 
composed  of  many  small  grains  or  fewer  larger  grains.  If  the  relationship  between  transport  rate  and 
grain  size  can  be  expressed  as: 

transport  rate  oc  D"  (7) 
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then,  for  the  ratios  in  Figure  1,  n  can  be  found  as: 


n  = 


log  (transport  ratio) 


(8) 


where  D\  is  the  smaller  grain  size  and  Di  is  the  larger.  Using  his  data,  King  (1991)  found  that  in  the 
sheetflow  regime:  n  ~  0.08,  as  shown  in  Figure  2,  demonstrating  the  small  influence  of  grain  size  on 
transport  rate. 


Figure  2.  Calculation  of  exponent  n  in  Equation  8  for  transport  ratios  shown  in  Figure  1 


In  the  IM  bed-load  regime,  whether  a  grain  in  the  surface  layer  moves  or  not  is  a  function  of  the  flow 
intensity  (t),  of  the  chance  of  impact  from  other  saltating  grains,  and  of  how  a  grain  is  sitting  and 
resting  against  its  neighbors.  At  sand  sizes  and  larger,  large  grains  are  less  easily  dislodged  and 
transported  than  small,  making  transport  a  function  of  grain  size. 

Sawamoto  and  Yamashita  (1986,  1987)  developed  a  relationship  for  sediment  transport  in  the  IM 
bed-load  regime.  Ignoring  the  effect  of  grain  size  on  the  friction  factor  (fw),  their  relationship 
indicates: 


_9 

transport  rate  oc  D 


(9) 
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The  description  of  the  Sleath  (1978)  experiments  indicate  that  both  grain  sizes  were  clearly  within 
the  IM  regime,  as  he  photographed  and  measured  the  motion  of  individual  particles.  As  seen  in 
Figure  2,  his  limited  data  set  also  indicates  that  for  this  regime,  n  should  be  near  -2.  There  is 
substantial  scatter  in  the  data  of  Kalkanis  (1964)  and,  unfortunately,  his  methodology  is  not  well 
described. 

The  threshold  of  motion  is  included  within  the  IM  regime  at  a  critical  Shields  parameter  value  of  the 
order  of  9C  ~  0.04  to  0.06.  Thus,  the  threshold  of  grain  motion  should  be  a  function  of  grain 
diameter.  Sleath  (1984)  (section  6.6. 1)  lists  many  common  threshold  of  motion  formulas  and  grain 
size  is  a  parameter  included  in  most  of  these  equations.  However,  most  of  these  formulas  indicate 
that  the  appropriate  grain  size  exponent  in  Equation  9  should  be  between  0  and  -1.  Clearly,  more 
data  would  be  useful  in  exploring  the  behavior  of  transport  in  this  regime  and  its  relationship  to  grain 
size. 

In  the  laboratory,  if  experiments  are  conducted  using  periodic  waves  to  transport  sand  at  rates  a  little 
above  threshold,  ripples  fonn  on  the  bed.  However,  these  ripples  take  time  to  form  (many  wave 
periods,  depending  upon  the  flow  intensity).  Prior  to  their  formation,  transport  is  in  the  IM  bed-load 
regime,  but  after  they  form,  it  is  in  the  ripple  regime.  The  situation  is  somewhat  more  complicated 
with  irregular  waves,  in  that  large  waves  will  occasionally  wipe  out  the  ripples,  causing  the  process 
to  start  over.  So  transport  will  alternate  between  these  regimes;  however,  it  will  usually  be  in  the 
ripple  regime.  As  shear  stresses  increase,  ripples  go  through  a  progression  of  modifications  until  they 
are  wiped  out  at  a  Shield’s  parameter  value  of  the  order  of  9  ~  0.8.  This  is  well  within  the  sheetflow 
bed-load  regime.  Thus,  the  presence  of  ripples  has  tended  to  mask  the  existence  of  two  bed-load 
regimes.  Studies  of  the  IM  bed-load  regime  must  generally  be  done  before  ripples  develop  or  at 
large-grain  diameters,  for  ripples  do  not  tend  to  form  for  grain  sizes  above  D  ~  0.8  mm. 

From  shallow-water  linear  wave  theory,  the  maximum  horizontal  velocity  (um)  is: 

^=\HWd  (io) 


where  H  is  wave  height  and  d  is  water  depth.  Using  this  relationship  plus  Equation  2,  the  boundary 
equation,  0  ~  0.2,  between  the  two  bed-load  regimes  is  found  to  be  equivalent  to: 


(p~pJ 

yfw 


(ii) 


where  y  is  the  breaking  criteria  (equal  to  Hbld).  Inserting  the  approximate  values  y  =  0.78,  /,',.  ~  0.015, 
and  ps  =  2.65  (quartz  sand  density),  yields  the  equivalent  boundary  relationship  in  a  breaking  wave: 
Hb!D  ~  230.  For  the  fine  sands  common  on  most  beaches  where  D  is  on  the  order  of  0.2  mm, 
breaking  wave  heights  are  almost  always  greater  than  a  few  centimeters;  so  this  value  is  almost 
always  exceeded,  and  bed-load  occurs  in  the  sheetflow  bed-load  regime  as  a  dense  layer  of  moving 
material,  many  grain  diameters  thick.  However,  on  coarse  grain  beaches  (gravel  or  shingle)  where  D 


8 


ERDC/CHL  CHETN-ll-48 
November  2005 

is  on  the  order  of  0.01  to  0.1  m,  this  Hb/D  value  is  rarely  exceeded  and  bed-load  occurs  in  the  IM 
regime  in  which  transport  decreases  as  grain  size  increases. 

Are  Two  Bed-load  Regimes  Necessary?  As  fine  grain  beaches  are  predominant  worldwide, 
and,  on  these,  the  presence  of  ripples  will  generally  render  calculations  of  transport  in  the  IM  bed¬ 
load  regime  unnecessary,  it  is  appropriate  to  consider  if  the  recognition  of  two  bed-load  regimes  can 
be  expected  to  lead  to  improved  estimates  of  sediment  transport  rates.  The  impact  is  likely  to  be  of 
significance  in  two  ways.  First,  on  coarse-grained  (gravel  and  cobble)  beaches,  suspended  load 
transport  does  not  occur  (the  suspension  parameter,  Equation  3,  is  large),  sheetflow  bed-load 
transport  does  not  occur  because  the  waves  are  not  energetic  enough,  and  ripples  do  not  generally 
form  at  grain  sizes  above  D  ~  0.8  mm.  Therefore,  by  default,  on  these  beaches  IM  bed- load  transport 
can  be  expected  to  be  the  dominant  transport  mechanism.  Secondly,  the  inclusion  of  two  regimes 
will  lead  to  improved  equations  for  bed-load  transport  in  the  sheetflow  regime.  The  description  of 
sheetflow  transport  should  be  relatively  simple.  It  is  approximately  independent  of  grain  size  and  a 
threshold  term  does  not  need  to  be  included,  because  threshold  of  motion  occurs  outside  of  this 
regime. 

It  is  not  mandatory  that  bed-load  transport  be  divided  into  two  regimes.  A  single  appropriate 
equation  could  be  devised  that  included  the  threshold  of  motion  and  the  behavior  described  in 
Figure  1 .  However,  the  types  of  formulations  commonly  employed  to  include  threshold  effects  have 
significant  limitations  (King  and  Seymour  1989).  For  example,  Sleath  (1978)  proposed  the  following 
transport  relationship: 

-iL  =  47(0-0,)^  (12) 

c oD 

where  Q  is  the  volume  transport  rate  per  unit  width  and  co  is  the  wave  frequency.  In  this  equation, 
the  threshold  tenn,  0C  (~  0.05)  substantially  affects  (depresses)  the  transport  rate  at  Shield’s 
parameter  values  well  above  0  ~  0.2.  At  0  =  0.7  the  transport  rate  is  only  90  percent  of  the  value  it 
would  be  if  the  threshold  term  were  not  included  in  the  equation.  Most  transport  equations  that 
include  a  threshold  term  include  it  in  a  similar  fashion. 

The  transport  equation  of  Madsen  and  Grant  (1976)  does  a  much  better  job  of  limiting  the  effect  of 
the  threshold  term  to  the  region  near  threshold: 

-^-  =  ce3  (13) 

wsD 

where  C  is  a  parameter  related  to  the  percentage  of  time  in  a  wave  cycle  when  the  velocities  are 
above  threshold.  The  C  values  are  defined  in  a  table  in  their  paper.  This  parameter  incorporates  the 
threshold  of  motion  because  C  =  0  at  0  =  0C;  it  reaches  its  full  value  of  C  =  12.5  at0  =  40c  ~0.2. 
However,  this  equation  does  not  treat  the  grain  size  dependence  correctly.  If fw  is  assumed  to  be 
independent  of  D,  and  ws  oc  D,  this  equation  implies  that  the  transport  is  proportional  to  D'1  for  all 
values  of  the  Shield’s  parameter. 
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Suspended  Load  Transport  Regime.  For  sediment  suspension  to  occur,  the  vertical  turbu¬ 
lence  must  be  sufficient  to  counteract  the  fall  velocity,  ws.  See  Nielsen  (1992)  for  a  detailed 
description  of  the  process.  The  level  of  turbulence  is  a  function  of  the  maximum  wave  induced 
bottom  shear  stress,  zm,  or,  alternatively,  the  friction  velocity,  it*.  This  reasoning  leads  directly  to  the 
Rouse  or  suspension  parameter,  Equation  3,  being  important  in  the  description  of  sediment 
suspension. 

Many  formulas  have  been  proposed  to  describe  suspended  load  sediment  transport.  Essentially  all  of 
these  indicate  that  the  rate  decreases  with  grain  size,  or  equivalently,  with  fall  velocity,  and  most 
indicate  that  the  transport  rate  is  strongly  inversely  dependent  on  grain  size.  For  a  typical  develop¬ 
ment,  the  suspended  load  sediment  transport  rate  within  a  vertical  column  of  unit  area  in  the  surf 
zone  can  be  given  as: 

r  d 

Q(x,t)  =  J  c(x,z,t )  x  u(x,z,t)dz  (14) 

where  c(x,z,t )  is  the  sediment  concentration  and  u(x,z,t )  is  the  sediment  velocity,  with  both  terms 
being  functions  of  cross-shore  (x)  and  vertical  (z)  position  and  time  ( t ). 

To  obtain  the  total  longshore  suspended  load  transport  rate,  Qys,  this  equation  can  be  integrated  over 
a  wave  cycle  and  across  the  surf  zone.  In  the  gradient  diffusion  model,  the  settling  of  grains  toward 
the  bed  at  the  fall  velocity,  ws,  is  balanced  by  turbulent  diffusion,  given  as: 

wsC(z)  =  -&^~  (15) 

dz 

where  s  is  the  eddy  diffusivity  of  sediment.  For  rippled  beds,  the  eddy  diffusivity  can  be  assumed 
constant  with  height  and  the  corresponding  concentration  profile  is  given  as: 

C(z)  =  C0e~z/I  (16) 

where  Co  is  the  reference  concentration  and  /  is  the  decay  length  scale. 

Numerous  equations  have  been  proposed  for  the  reference  concentration,  the  decay  length  scale,  and 
for  the  form  of  the  concentration  profile.  Nielsen  ( 1 992)  proposed  that  the  reference  concentration  is 
proportional  to  the  cube  of  the  ripple  Shield’s  parameter,  and  the  decay  length  is  a  function  of  ripple 
geometry.  This  development  leads  to  the  suspended  load  transport  rate  being  approximately  equal  to 
D'3.  Transport  relationships  developed  for  flat  beds  generally  give  similar  results. 

LONGSHORE  SEDIMENT  TRANSPORT  MODELS:  FourLST  models  are  discussed  later  to 
show  how  they  incorporate  the  physics  of  grain-size  variation  and  its  effect  on  the  transport  rate:  the 
CERC  formula,  the  Kamphuis  model,  the  Bailard  model,  and  the  Madsen  and  Tajima  model.  The 
models  selected  are  not  intended  to  be  inclusive,  but  merely  representative  of  classes  of  models.  The 
CERC  equation,  containing  one  term  for  the  calculation  of  total  load  (combined  bed-load  and 
suspended  load),  is  the  simplest  formula  in  general  use.  The  Kamphuis  formula  is  also  a  one -term, 
total  load  model,  but  explicitly  includes  the  effects  of  wave  period,  beach  slope,  and  grain  size.  The 
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Bailard  formula  is  representative  of  models  that  partition  the  transport  into  bed-load  and  suspended 
load  transport.  The  Madsen  and  Tajima  model  is  representative  of  the  complex  computer  routines 
that  model  the  hydrodynamics  and  sediment  dynamics  stepwise  across  the  surf  zone,  and  thus  predict 
not  only  the  total  bed-load  and  suspended  load  LST,  but  also  its  cross-shore  distribution. 

CERC  Formula.  The  most  common  method  of  calculating  the  total  LST  rate  in  use  today  is  the 
CERC  formula  {Shore  Protection  Manual  1977, 1984;  Coastal  Engineering  Manual  2002),  which  is 
given  as: 


/  = 


K 

!6Vy 1 


3  5 

P g2Hl  sin(2a/j) 


(17) 


where  /  is  the  total  immersed  weight  LST  rate,  K  is  an  empirical  coefficient,  Hb  is  the  breaking  wave 
height,  jh  is  the  breaker  index  (typically  equal  to  0.78),  and  a./,  is  wave  angle  at  breaking.  An  equiva¬ 
lent  expression  is: 


0  = 


K 


16 


ix~p) 


-Hi  sin(2a/;) 
Y 


(18) 


where  Q  is  the  total  bulk  volume  LST  rate  (volume  of  the  grains  plus  voids). 

This  model,  which  is  based  on  the  assumption  that  the  total  longshore  sediment  transport  rate  is 
proportional  to  longshore  energy  flux,  was  developed  from  the  pioneering  work  of  Bagno  Id  in  the 
early  1960s  and  further  developed  by  Komar  and  Inman  (1970).  Lor  a  detailed  description  of  the 
development,  see  Dean  and  Dalrymple  (2002). 

The  limitations  of  the  CERC  formula  are  well  known.  With  adequate  calibration,  the  CERC  formula 
can  estimate  LST  rates  within  +  50  percent.  However,  without  calibration,  the  CERC  fonnula  only 
provides  an  accuracy  of  one  to  two  orders-of-magnitude  (Greer  and  Madsen  1978;  Lowler  et  al. 
1995;  Wang  et  al.  1998).  However,  the  CERC  equation  continues  to  be  useful,  primarily  because  of 
its  simplicity  and  because  of  the  failure  of  more  sophisticated  models  to  clearly  demonstrate  substan¬ 
tially  superior  accuracy  relative  to  the  effort  required  to  employ  them. 

Variation  of  CERC  K  With  Grain  Size.  The  recommended  value  of  the  CERC  K  coefficient 
depends  upon  whether  the  breaking  wave  height  is  specified  in  terms  of  the  significant  wave  height 
or  the  root-mean-square  wave  height.  See  Galvin  and  Scheppe  (1980)  for  a  discussion  of  the 
relationship  between  these  quantities.  Following  Schoonees  and  Theron  ( 1 994),  in  this  technical  note 
K  will  be  subscripted,  where  necessary,  to  reduce  ambiguity.  K{sig)  will  refer  to  K  values  associated 
with  the  significant  wave  height,  and  K(rms)  will  refer  to  root-mean-square  wave  height  K  values. 
Where  necessary  to  convert  between  these  quantities,  the  equation:  K(sig)  =  Vi  Kfrms)  is  used.  The 
Shore  Protection  Manual  (1984)  recommends  a  value  of  K(sig)  =  0.39  (K(rms)  =  0.78),  which  was 
based  upon  41  measurements  of  LST  rate. 
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Several  attempts  have  been  made  to  modify  the  K  coefficient  to  improve  the  CERC  formula  and  this 
has  frequently  involved  relating  K  to  the  grain  size,  as  most  investigators  have  concluded  that  the 
LST  rate  should  decrease  as  grain  size  increases  (Dean  1978;  Walton  and  Chui  1979;  Bruno  et  al. 
1980;  Dean  et  al.  1982;  Dean  1989;  Dean  and  Dalrymple  2002). 

Komar  (1988)  examined  most  of  the  available  field  data,  and  questioned  the  validity  of  some.  He 
concluded  that  existing  data,  most  of  which  were  for  sand-size  beaches,  showed  little  relationship 
between  K  and  grain  size.  Theory  and  data  showed  a  trend  of  decreasing  K  with  increasing  grain 
size,  but  the  data  were  too  limited  to  establish  a  correlation.  Komar  stressed  that  K  should  depend  on 
sediment  grain  size,  and  the  absence  of  such  a  trend  in  his  analysis  resulted  from  the  imperfect 
quality  of  the  data.  For  sand-sized  grains,  he  recommended  a  K(rms)  =  0.57,  based  upon  70  data 
points.  In  comparison,  Kraus  et  al.  (1982)  used  25  data  points  to  recommend  a  value  of  K(rms)  =  0.58. 
Bodge  and  Kraus  (1991)  gave  another  extensive  review  of  problems  associated  with  use  of  the 
CERC  equation. 

Schoonees  and  Theron  (1993)  compiled  and  reviewed  35  data  sets  containing  273  measurements  of 
bulk  longshore  transport  rate.  They  devised  a  rating  system  and  assessed  the  quality  of  the  sets.  Then 
Schoonees  and  Theron  (1994)  used  portions  of  the  data  set  to  evaluate  several  version  of  the  CERC 
formula.  They  also  divided  the  data  into  those  with  grain  sizes  less  than  1  mm  and  those  with 
greater.  For  the  fine-grained  data  set  (206  points),  they  found  a  best  fit  K(sig)  =  0.20.  If  only  the 
highest  quality  data  sets  were  used  (as  per  Schoonees  and  Theron  1 993 , 46  data  points)  they  obtained 
a  best  fit  K(Sig)  =  0.41.  However,  the  correlation  was  only  slightly  higher  for  the  second  set  than  the 
first  ( R2=  coefficient  of  determination  =  0.77  for  second  set  versus  0.72  for  the  first).  For  the  data  set 
with  sediment  larger  than  1  mm  (34  data  points),  they  obtained  K(sig)  =  0.01,  but  the  correlation 
coefficient  was  only  i?2  =  0.11. 

It  is  generally  accepted  that  there  is  a  relationship  between  the  CERC  K  term  and  grain  size. 
Formulas  relating  those  two  quantities  have  been  proposed  by  Swart  (1976),  Bailard  (1984), 
USAEWES  (1985),  Kamphuis  et  al.  (1986)  and  del  Valle  et  al.  (1993). 

Swart  (1976)  used  a  different  form  of  the  CERC  relationship.  Following  Schoonees  and  Theron 
(1994),  his  relationship  can  be  expressed  as: 

V>)  =  U51°s4°'001%5J  <19) 


where  the  median  grain  diameter  is  in  meters.  The  relationship  is  valid  between 
0. 1  mm  <  D50  <1.0  mm.  He  commented  that,  intuitively,  he  would  have  expected  a  steeper  curve. 
The  Bailard  (1984)  relationship  is  given  below  in  the  development  of  his  model.  USAEWES  (1985) 
gives  a  relationship  that  relates  K  to  the  fall  velocity  and  also  to  the  wave  height.  This  expression, 
which  was  calibrated  with  the  field  data  available  at  that  time,  is: 


K(sig)  =  0.1637  log 


f  rr  \ 

gHb 


\Ws  J 


-0.0773 


(20) 
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Kamphuis  et  al.  (1986)  developed  a  relationship  based  upon  dimensional  analysis,  which  includes 
the  wave  height  and  breaker  index,  as  well  as  the  grain  size.  His  expression  is: 


h*  )  =  M22^Hb/D  <21) 

del  Valle  et  al.  (1993)  developed  an  empirically  based  relationship  for  K,  which  shows  decreasing 
values  of  K  with  larger  grain  sizes.  The  equation  was  based  on  data  presented  by  Komar  (1988)  and 
data  obtained  from  their  research  along  the  Adra  River  Delta,  Spain.  The  relationship  is: 

Kj  ,  =  1.4e('2'5Ao)  (22) 

yrms)  v  7 


where  D50  is  expressed  in  millimeters.  This  relationship  is  strongly  dependent  on  the  relatively  larger 
median  sand  grain  sizes  from  the  Adra  River  Delta  (D50  =  0.44  to  1.5  mm),  and  they  cautioned 
against  its  use  for  grain  sizes  above  1.5  mm. 

Figure  3  shows  data  relating  the  CERC  K  term  to  sediment  grain  size.  It  includes  the  field  studies 
listed  in  Komar  (1988),  plus  the  gravel  and  shingle  transport  studies  ofNicholls  and  Webber  (1987), 
Chadwick  (1989),  andNicholls  and  Wright  (1991)  from  England,  and  Leidersdorf  et  al.  (1988)  from 
Alaska.  Though  there  is  a  clear  trend  of  decreasing  transport  with  increasing  grain  size,  there  is  also 
extensive  scatter  in  the  data,  with  more  scatter  for  the  larger  grain  sizes. 

The  trend  in  the  data  suggests  a  relationship  of  the  form: 


_  0.1  mm 

K(«g)  ~  — 5~ 


(23) 


However,  due  to  the  wide  divergence  of  the  data,  this  formula  should  only  be  used  for  initial 
estimates,  and  should  not  be  substituted  for  site-specific  field  measurements.  The  fonnulas  of  Swart 
(1976)  and  del  Valle  et  al.  (1993)  are  also  shown  over  the  ranges  for  which  the  authors  claim 
validity.  For  sand  size  materials,  all  three  equations  predict  similar  K  values.  However,  for  larger 
gravels,  their  fonnulas  predict  a  rapidly  decreasing  K,  which  is  outside  the  scatter  of  the  data,  and 
thus,  their  fonnulas  should  only  be  employed  within  their  stated  data  ranges. 


Kamphuis  Model.  Kamphuis  (1991)  expanded  upon  his  earlier  work  leading  to  Equation  20  and 
developed  a  relationship  for  estimating  longshore  sediment  transport  rates  based  upon  dimensional 
analysis  and  calibrated  this  relationship  using  physical  model  experiments.  The  equation,  which 
Kamphuis  (2002)  found  to  be  applicable  to  both  field  and  model  data,  is  expressed  as: 

S,  =  2.nH]bTYmfsDf»  sin“(2 a„ )  (24) 

in  which  Qu  is  the  transport  rate  of  underwater  mass  (kg/sec),  Tp  is  the  peak  wave  period  (sec),  /«/,  is 
the  beach  slope  from  the  breaker  line  to  the  shoreline,  and  D50  is  the  median  grain  size  (mm). 
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D50,  median  grain  size  (mm) 


Figure  3.  Relationship  between  CERC  K  term  and  grain  size 

Kamphuis  (2002)  uses  the  same  equation,  but  redefines  the  beach  slope  as  the  slope  that  causes 
breaking,  i.e.,  the  slope  over  one  or  two  wavelengths  offshore  of  the  breaker  line.  Equation  24  is 
appealing  because  it  includes  the  effects  of  wave  period,  which  influences  wave  breaking  (Galvin 
1968),  beach  slope,  and  grain  size.  Numerous  other  researchers  have  developed  formulas  that  relate 
the  LST  rate  to  a  suite  of  parameters  that  include  grain  size.  Some  of  these  are  van  Hijum  and 
Pilarczyk  ( 1 982);  Brampton  and  Motyka  ( 1 984);  Morfett  (1988);  Chadwick  (1989);  and  Van  Wellen 
et  al.  (2000). 

Bailard  Model.  Bailard  (1981)  developed  an  energy  based  surf  zone  sediment  transport  model 
based  upon  Bagnold’s  (1963,  1966)  steady  flow  models.  His  model  had  both  bed-load  and  sus¬ 
pended  load  components.  Both  components  were  expressed  in  terms  of  various  instantaneous 
velocity  components,  which  limited  its  usefulness.  Bailard  (1984)  made  several  simplifying 
assumptions  and  found  that  he  could  collapse  both  the  bed-load  and  suspended  load  components  to 
the  form  of  the  CERC  equation  if  the  K  term  was  of  the  form: 

K,  ,  =  0.05  +  2. 6 sin2  ( 2^ )  +  0.0096m  +  0.007^-  (25) 

1  ' 1  ws 

where  ws  is  the  fall  speed  of  the  sediment,  m  is  the  beach  slope,  and  is  the  maximum  oscillatory 
velocity  magnitude  at  breaking  obtained  from  Equation  10.  In  Equation  25,  the  first  three  terms  are 
contributions  from  the  bed-load  portion  of  his  model,  and  the  last  is  the  suspended  load  portion. 
Also,  the  m  term  is  usually  neglected  because  of  its  small  size. 
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Madsen  and  Tajima  Model.  Madsen  and  Taj ima  (Taj ima  2004)  developed  a  computer  routine  to 
model  surf  zone  sediment  transport.  The  code  is  in  the  form  of  two  programs  that  run  sequentially:  a 
hydrodynamics  model  and  a  sediment  transport  model.  The  hydrodynamic  model  calculates  the 
forcing  functions  needed  to  drive  the  sediment  transport  model  at  each  point  on  the  profile,  and 
includes  modules  for  nonlinear  wave  propagation,  wave  breaking,  surface  rollers,  and  nearshore 
currents.  The  sediment  transport  model  calculates  the  transport  at  each  profile  point  and  includes 
bed-load  and  suspended  load  modules. 

Model  inputs  include  wave  height,  period,  and  direction,  median  sediment  size  or  size  distribution, 
beach  profile,  and  various  physical  and  model  control  parameters.  The  model  contains  no  adjustable 
calibration  coefficients  such  as  the  CERC  K  coefficient.  Model  outputs  include  wave-average  values 
for  wave  height,  setup,  and  current  across  the  profde,  the  cross-shore  distribution  of  bed,  suspended, 
and  total  load  LST  rate,  along  with  bulk  values  for  each,  and  vertical  sediment  distributions  at  user 
defined  cross-shore  locations. 

An  example  of  the  suspended  load  and  bed-load  model  predictions  is  shown  in  Figure  4.  Q  is  the 
bulk  volume  (sediment  plus  voids  (p  =  0.4))  yearly  LST  rate  that  would  occur  for  a  constant 
breaking  wave  height  with  a  period  of  (T  =  8  sec)  and  breaking  wave  angle  of  (a  =  15  deg)  on  a 
grain-size-adjusted  equilibrium  beach  profde.  The  suspended  load  transport  shows  a  significant 
decrease  in  transport  rate  with  increasing  grain  size,  while  the  bed-load  rate  is  largely  grain  size 
independent. 


Grain 

Sizes 

(mm) 


-0.13 

— 0- 

-0.18 

A 

-0.23 

— B- 

-0.28 

-0.33 

Significant  Wave  Height  (m) 


Figure  4.  Madsen  and  Tajima  model  suspended  load  (solid  lines)  and  bed-load  (dashed  lines)  LST  rate 
predictions 
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MODEL  LST  PREDICTIONS  OF  EFFECTS  OF  GRAIN  SIZE:  To  show  how  variations  in 
grain  size  affect  the  model  predictions,  the  four  models  previously  described  were  run  using  two 
typical  beach  sediments  with  grain  sizes  D  =  0.13  and  0.33  mm.  The  models  used  wave  periods  of 
T=  8  sec,  breaker  wave  angles  of  op>  =  15  deg,  and  equilibrium-shaped  profiles  for  each  grain  size. 
The  results  for  each  model  are  shown  in  Figure  4  as  ratios  of  the  LST  rate.  The  ratios  plotted  are  the 
smaller  grain  size  LST  rate  (D  =  0.13  mm)  divided  by  the  larger  (D  =  0.33  mm)  as  a  function  of 
significant  wave  height. 

It  was  previously  shown  (Equation  1 1)  that  the  ratio  Hh/D  is  approximately  equivalent  to  the  Shield’s 
parameter,  9.  As  the  grain  sizes  in  Figure  5  are  held  constant,  the  horizontal  axis  on  this  figure  is 
analogous  to  the  horizontal  axes  in  Figures  1  and  2.  Both  9  and  Hb/D  can  be  thought  of  as  ratios  of 
the  tangential  force  mobilizing  the  sediment  to  the  gravitational  force  resisting  the  motion.  In 
Figure  5,  all  bed-load  transport  should  be  considered  to  be  in  the  sheetflow  regime,  because  the 
boundary  for  the  larger  grain  size  should  be  at  a  wave  height  of  approximately  7  cm.  (None  of  the 
models  explicitly  include  a  term  representing  IM  bed-load  transport.) 


Significant  Wave  Ht  (m) 


Figure  5.  Comparison  of  model  predictions  for  ratio  of  0.13  mm  grain  size  to  0.33  mm  transport  rate 

The  CERC  formula  and  the  Bailard  bed-load  formula  have  no  grain  size  dependence,  so  for  both,  the 
ratio  =  1 .  The  Kamphuis  total  load  ratio  and  the  Madsen  and  Tajima  bed-load  ratio  are  both  slightly 
less  than  one.  For  the  Madsen  and  Tajima  model,  this  is  due  to  the  larger  bed  roughness  for  the 
larger  grain  size.  The  Kamphuis  ratio  is  not  constant  because  equilibrium  beach  slopes  were  used  for 
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each  grain  size,  and  the  Kamphuis  formula  explicitly  includes  beach  slope.  If  the  same  profile  had 
been  used  for  both  Kamphuis  runs,  the  transport  ratio  would  be  1 .26,  implying  greater  transport  for 
the  smaller  grain  size.  The  ratio  for  the  Bailard  suspended  load  model  is  the  inverse  of  the  fall 
velocity  ratio  for  the  two  grain  sizes  (=  3.29).  The  Madsen  and  Tajima  suspended  load  model 
includes  a  fonnulation  of  the  type  discussed  in  the  section  on  the  Suspended  Load  Regime  above, 
and,  as  expected,  shows  significantly  greater  transport  for  the  smaller  grain  size. 


Table  2  shows  the  value  of  the  exponent  n  for 
each  model  using  the  simple  assumption  that 
transport  is  proportional  to  Dn  (Equation  7). 
These  n  values  were  either  taken  directly  from 
the  model  formulation  or  were  calculated  using 
Equation  8.  In  general,  the  table  shows  that 
bed-load  transport  is  largely  independent  of 
grain  size,  while  suspended  load  transport 
strongly  decreases  as  grain  size  increases. 


CERC  FORMULA  PARADOX:  Figure  6  shows  the  ratio  of  suspended  load  transport  to  bed-load 
transport  as  a  function  of  wave  height  and  grain  size  as  predicted  by  the  Madsen  and  Tajima  model 
(wave  period  T=  8  sec,  breaker  wave  angle  otb  =  15  deg,  and  equilibrium-shaped  profile  for  each 
grain  size).  This  model  predicts  what  most  researchers  feel  is  true,  that  suspended  load  transport  is 
usually  much  more  significant  than  bed-load  transport  on  most  fine-grain  beaches  most  of  the  time. 
The  discussion  indicates  that  suspended  load  transport  has  a  strong  dependence  on  grain  size. 
However,  this  would  indicate  that  the  CERC  formula,  which  has  no  dependence  on  grain  size,  or  the 
Kamphuis  formula,  which  has  only  a  weak  dependence,  should  not  have  much  predictive  skill.  This 
contradiction  has  been  termed  “the  CERC  formula  paradox”  by  Nielsen  (1988,  1992). 

This  paradox  is  at  least  partially  resolved  by  noting  that  grain  size  is  related  to  beach  slope,  as  seen 
in  the  equilibrium  beach  profile  concept  (Dean  1977,  1991),  and  LST  also  is  a  function  of  beach 
slope,  m,  and  the  surf  similarity  parameter,  2: 


Table  2 

Model  Dependence  on  Grain  Size 

Model 

Transport  =  D" 

CERC  Total  Load 

D° 

Kamphuis  Total  Load 

D-0.25 

Bailard  Bed-load 

D° 

Bailard  Suspended  Load 

w;' 

M  &  T  Bed-load 

Dn  where  0.2<  n  <0.4 

M  &  T  Suspended  Load 

D"  where  -3.6<  n  <2.5 

5  = 


(26) 


Larger  grain  sizes  (shown  previously  to  reduce  the  transport  rate)  are  associated  with  beach  profiles 
having  steeper  slopes,  but  steeper  slopes  are  associated  with  greater  transport  rates.  David  Altman 
and  Robert  Dean  have  looked  at  the  LST  rate  associated  with  beach  fills  and  found  greater  transport 
on  initial  steep  beach-fill  profiles  than  on  later,  less  steep  profiles.  (Personal  Communication,  20 
June  2005,  Robert  Dean,  Professor,  University  of  Florida).  The  surf  similarity  parameter  has  been 
related  to  breaker  type  (Battjes  1974),  with  milder  slopes  producing  spilling  breakers  and  steeper 
slopes  producing  plunging  breakers.  Nielsen  (1991)  commented  that  the  type  of  breaking  influences 
the  surf  zone  turbulence  levels,  with  mild  slope  spilling  breakers  producing  less  near-bed  turbulence, 
and  thus  decreasing  the  suspended  transport.  The  beach  slope  tenn  is  explicitly  included  in  the 
Kamphuis  model,  with  the  effects  on  transport  discussed  previously.  Several  other  researchers 
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Grain 

Sizes 

(mm) 


-m- 

0.13 

-e- 

0.18 

0.23 

□ 

0.28 

0.33 

Significant  Wave  Height  (m) 


Figure  6.  Ratio  of  suspended  load  to  bed-load  LST  rate  as  predicted  by  Madsen  and  Tajima  model 

(Kamphuis  and  Readshaw  1970;  Kamphuis  and  Sayao  1982;  Ozhan  1982)  have  discussed  relation¬ 
ships  between  the  LST  rate  and  the  surf  similarity  parameter.  This  relationship  is  further  explored  by 
Smith  et  al.  (2004). 

The  implication  is  that  the  beach  slope,  which  is  partially  a  function  of  grain  size,  should  tend  to 
counteract  the  decreased  LST  rate  associated  with  the  increased  grain  size  effects  discussed  above. 
However,  it  is  clear  from  the  trend  shown  in  Figure  3  that  the  dominant  effect  of  increased  grain  size 
is  to  decrease  the  LST  rate. 

SUMMARY  AND  CONCLUSIONS:  What  is  classically  considered  the  bed-load  transport  regime 
is  actually  composed  of  two  regimes,  an  Initial  Motion  (IM)  bed-load  regime,  and  a  sheetflow  bed¬ 
load  regime. 

For  coarse  cobble  and  shingle  beaches,  where  the  top  layer  of  the  bed  is  not  fully  mobilized  by  wave 
action,  the  dominant  regime  is  the  IM  bed-load  regime.  Though  not  well  studied,  the  transport  rate  is 
inversely  related  to  grain  size,  with  n,  the  grain  size  exponent  in  Equation  7,  of  the  order  of  -0.5 
to  -2. 


In  the  sheetflow  bed-load  regime,  transport  is  largely  independent  of  grain  size. 

In  the  suspended  load  regime,  transport  is  strongly  inversely  related  to  grain  size  with  n  being  on  the 
order  of  -3  (without  considering  the  effects  of  beach  slope). 
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In  the  ripple  regime,  transport  is  expected  to  show  the  same  grain  size  dependence  as  in  the 
suspended  load  regime.  In  the  absence  of  evidence  to  the  contrary,  in  the  swash  regime,  the 
dependence  of  the  longshore  transport  on  grain  size  is  expected  to  be  similar  to  that  in  bed-load  or 
suspended  load,  depending  upon  which  is  dominant. 

A  relationship,  Equation  23,  between  the  grain  size  and  the  K  coefficient  in  the  CERC  formula  is 
presented.  However,  there  is  substantial  scatter  in  the  data  and  the  reader  is  cautioned  to  calibrate  the 
value  of  K  using  measured  transport  rates  if  the  CERC  formula  is  to  be  employed. 

A  series  of  LST  rate  models  is  presented  and  their  predictions  for  the  variation  of  transport  with 
grain  size  are  presented  in  Table  2. 

The  assumption  that  the  dominant  mode  of  transport  on  most  beaches  is  in  the  suspended  load 
regime  leads  to  a  paradox.  In  general,  LST  rate  data  do  not  show  the  strong  inverse  dependence  on 
grain  size  that  would  be  expected  in  sediment  suspension.  This  paradox  is  at  least  partially  resolved 
by  noting  that  there  is  a  relationship  between  grain  size  and  beach  slope,  and  that  variations  in  beach 
slope  can  also  modify  the  transport  rate. 

In  general,  an  increase  in  the  median  grain  size  will  decrease  LST  rates  in  the  surf  zone.  If  a  simple 
exponential  relationship  between  transport  rate  and  grain  size  ( D )  is  needed  (Equation  7:  LST  rate  oc 
D"),  the  most  appropriate  value  for  the  exponent  n  should  be  of  the  order  of  -1,  as  seen  from 
Equation  23.  However,  this  technical  note  argues  that  this  is  clearly  a  simplistic  view  of  surf  zone 
sediment  dynamics.  A  more  realistic  (though  still  highly  simplified)  approach  would  be  that,  for  fine 
grain  sediments  (D50  on  the  order  of  0. 15  to  0.30  mm),  suspended  load  transport  should  dominate 
and  n  should  be  somewhere  within  the  range  of  -0.5  to  -3.0.  For  coarse  sands  (D50  around  1 .0  mm), 
sheetflow  bed-load  transport  should  dominate  (Figure  6)  and  the  transport  rate  should  be  nearly 
independent  of  grain  size  (n  =  0).  For  large  gravel  and  shingle  (D$o  >  20  mm),  the  dominant  transport 
should  be  in  the  IM  bed-load  regime,  with  n  within  the  range  of  -0.5  to  -2.0.  To  state  it  another  way, 
the  exponent  n  in  Equation  7  is  itself  a  function  of  grain  size. 

ADDITIONAL  INFORMATION:  Questions  about  this  CHETN  can  be  addressed  to  Dr.  David  B. 
King,  Jr.  (Voice:  601-634-2087),  (Fax:  601-634-4314),  (e-mail:  David.B.King@erdc.usace 
.army.mil),  U.S.  Army  Engineer  Research  and  Development  Center,  Coastal  and  Hydraulics 
Laboratory,  HF-CI,  3909  Halls  Ferry  Road,  Vicksburg,  MS  39180-6199.  The  contributions  of 
Ms.  Julie  Rosati,  Dr.  Bradley  Johnson,  and  Mr.  Ernest  Smith,  are  gratefully  acknowledged. 
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